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We show that dark solitons in ID Bose liquids may be created by absorption of a single quanta 
of an external ac field, in a close analogy with the Einstein's photoelectric effect. Similarly to the 
von Lenard's experiment with photoexcited electrons, the external field's photon energy huj should 
exceed a certain threshold. In our case the latter is given by the soliton energy e s {hq) with the 
momentum hq, where q is photon's wavenumber. We find the probability of soliton creation to 
have a power-law dependence on the frequency detuning u) — e s /h. This dependence is a signature 
of the quantum nature of the absorption process and the orthogonality catastrophe phenomenon 
associated with it. 
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I. INTRODUCTION 

The existence of dark solitons (DS) is among the most 
spectacular manifestations of the role played by weak 
inter-particle interactions in ID cold atomic gases 
Such solitons are macros comically large areas of partially, 
or even completely depleted gas, which propagate co- 
herently without any dispersion. It is natural to in- 
terpret these objects as localized solutions of the semi- 
classical Gross-Pitaevskii equation [2| . Correspondingly, 
the means to create DS, employed so far, required a 
macroscopic classical perturbation applied to the atomic 
cloud. An example of the latter is the phase imprinting 
technique [3[, where a finite fraction of the ID atomic 
cloud is subject to an external potential for a certain 
time. Once the potential is switched off and the gas is 
allowed to evolve, the DS is formed around the place with 
the maximal gradient of the potential. 

Drawing an analogy with the electronic field-emission 
from a metal: a pulse of a strong external electric field 
may lead to creation of free electrons outside of the metal 
surface. It is well-known, however, from the time of von 
Lenard and Einstein 0, [fl that this is not the only way 
to excite electrons. Indeed, a weak ac field results in a 
photoelectric current, as long as the energy of its quanta 
exceeds the threshold given by the work function of the 
metal. The difference between the field-emission and the 
photoelectric effects is that the latter essentially utilizes 
the quantum nature of the electromagnetic radiation. Is 
there an analog of the photoelectric effect for excitation of 
DS? Namely can the DS be created by a weak ac radiation 
with the frequency exceeding a certain threshold? 

In the framework of the Gross-Pitaevskii equation the 
answer on these questions is negative. Indeed, a weak ex- 
ternal field may lead to excitation of the linear waves, if 
its wavenumber and frequency satisfy Bogoliubov disper- 



sion relation, but not to creation of DS. However, treating 
the Bose liquid beyond the semiclassical Gross-Pitaevskii 
approximation reveals that creation of DS in response 
on an absorption of a single quanta with an above-the- 
threshold energy is actually possible. In analogy with the 
photoelectric effect we call this phenomenon the photo- 
solitonic effect. The threshold energy is given by the 
energy of DS £ s {p) with the momentum p = hq, where q 
is the photon wavenumber. Creation of DS requires an 
ac field with frequency Hu> > s s (p). Notice that no com- 
parison of the external frequency lo and the DS energy e s 
ever appears in the Gross-Pitaevskii treatment. 

Consider ID Bose liquid subject to an external ac po- 
tential with the wavenumber q and frequency lo. Such 
a radiation may be created using Bragg scattering tech- 
nique [a, 0]- In these experiments the ac potential has 
been created by the interference pattern of two non- 
collinear optical beams with the differential frequency 
lo and the x-component of the differential wavevector q, 
Q. At zero temperature, according to the fluctuation- 
dissipation theorem, the probability to absorb the radi- 
ation is given by the dynamic structure factor (DSF), 
defined as the density-density correlation function 

S(q,u) = fdxdte^-^ (p(x,t)p(p,0)) , (1) 

where p(x, t) is the density operator. Rewriting DSF in 
the Lehman representation in terms of exact many-body 
eigenstates of the system \n) with energies e„, one finds 



S(<l,u) = ^2\(n\p q \0)\ 2 5(e n - e - fuo) 



(2) 
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where p q is the Fourier component of the density and n = 
corresponds to the ground state. Since the momentum 
is a good quantum number, only the many-body states 
with the total momentum p = hq contribute to the sum 
in the r.h.s. of Eq. ©. 

For the model with the short-range repulsive interac- 
tions the many-body spectrum has been evaluated ex- 
actly using the Bethe ansatz (BA) method 0. Lieb 
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has identified two characteristic modes in the excitation 
spectrum of the model [9(, known as Lieb I and II modes 
with the dispersion relations £1,2 (p), Fig- HI The two are 
given correspondingly by the particle and hole excita- 
tions in the set of the BA quasi-momenta. The hole-like 
mode £1 (p) is shown to be the lower bound of the many- 
body spectrum with a given momentum p. According to 
Eq. @ absorption is only possible if fao > £2(^0) ■ Em- 

ting a numerical implementation of the algebraic BA 
Caux and Calabrese [ll[ have shown that DSF is 
indeed non-zero for all energies in excess of £2(faq) and 
is peaked at the particle- like mode fao = Si(hq). In the 
limit of the weakly interacting gas the latter approaches 
the Bogoliubov dispersion relation 0, [12] 



£i(hq) — > fao B (q) = vshq^/l + {hq/2mv B ) 2 , (3) 

where vb is the Bogoliubov sound velocity and m is the 
boson mass. 

This observation offers a way to interpret absorption 
in a vicinity of the Lieb I mode, £i(p), in terms of 
weakly interacting Bogoliubov quasiparticles. Consider, 
e.g., a photon with some wave vector q and energy fao 
slightly below the value £i(hq), i.e., fao < £i(hq), see 
point A in Fig. [TJa). The energy and momentum con- 
servation laws allow for such photon to create two Bo- 
goliubov quasiparticles, fao = £i(Hq — p) + £%{p)- In 
the limit £\{hq) — fao <C mv B , one of the two parti- 
cles has small momentum and may be viewed as a "soft" 
phonon. For smaller initial photon energies, the result- 
ing two quasiparticles split the photon momentum more 
evenly, until the photon energy reaches the limiting value 
fao = 2£i(hq/2). If the photon energy is decreased below 
this threshold, a creation of more than two quasiparticles 
is needed to satisfy the conservation laws. Upon further 
lowering fao, more quasiparticles are created in the pro- 
cess of photon absorption. Once the photon energy fao 
approaches the line e — vufaq, Fig.[TJ the energy and mo- 
mentum of the absorbed photon is split between infinitely 
many soft phonons. 

Below this line the described process of dividing the 
energy and momentum between the quasiparticles does 
not work any more. Nevertheless the many-body spec- 
trum persists down to the lower value £2(faq) < vsfaq 
and the algebraic BA calculations ll| show that there is 
a finite absorption probability in the energy window 



£2{faq) < fao < vsfaq. 



(4) 



What is the absorption mechanism in this window, where 
the conservation laws forbid excitation of any number of 
quasiparticles or phonons? 

The clue to answer this question appeared in the 1976 
paper of Kulish, Manakov and Faddeev [l3[ , who noticed 
that the hole-like Lieb II mode approaches dispersion re- 
lation of DS in the weakly interacting limit 
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FIG. 1: Momentum-energy plane for excitations. A photon 
is represented by a solid dot. (a) Absorption of a photon (dot 
marked by ^4) with energy and momentum slightly below the 
quasiparticle spectrum £i(p), leads predominantly to creation 
of two quasiparticles with momenta q — p and p (the latter 
is determined by the shown geometrical construction). As A 
approaches the line e — vbP, the number of excited quasi- 
particles increases, (b) Once photon energy-momentum (dot 
marked B) fall below the line e = VbP, its absorption involves 
creation of a soliton. Other excitations created in the course 
of absorption may be treated as phonons if B is close to the 
boundary £2(f>). The shown "sound cone" with Bogoliubov 
velocity vb determines the range of possible momenta of a 
generated dark soliton, p- < p s < p+. 



It means that the many-body states with the energy in 
the vicinity of £2 must be viewed as quantized DS parti- 
cles. Correspondingly the photon absorption in the en- 
ergy window ^ necessarily involves excitation of DS 
along with Bogoliubov quasiparticles and/or phonons. 
Consider e.g. a photon with the energy immediately 
above the Lieb II mode, £2 < fao < £2 + mv 2 B (point 
B in Fig. Q] (b)). Drawing the "sound cone" with the 
slope vb down to the intersections with £2(p), one finds 
the range of the possible momenta of DS 



p- < p s < p+ 



(0) 



Slip) 



(5) 



which satisfy the conservation laws. Indeed, DS with the 
momentum p_ accompanied by a phonon, propagating 
in the direction of the external momentum hq, obviously 
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satisfies the energy and momentum conservation. Simi- 
larly, DS with the momentum p + must be accompanied 
by the counter-propagating phonon. Any other soliton 
from the momentum window © requires excitation of a 
certain superposition of the forward and backward prop- 
agating phonons. 

In this paper we evaluate probability W quj (p s ) to ex- 
cite DS with the momentum p s in the range ^ upon 
absorption of a photon with the wavenumber q and fre- 
quency uj > S2{f>q)/h. We show that such a probability 
is heavily shifted towards the lower boundary of the in- 
terval p s — p^(q,uj), i.e. DS is preferentially excited 
along with the forward moving phonon. At larger pho- 
ton energies, while still in the interval ([4]), DS is excited 
with highest probability along with the forward mov- 
ing Bogoliubov quasiparticle. Its energy and momen- 
tum may be found geometrically by plotting the replica 
of the Bogoliubov dispersion curve which starts at some 
point along the Licb II mode, £2(^)7 and passes through 
the point (hq, huj) representing external photon. We also 
show that the total probability to excite any DS scales as 
a power of the blue detuning from the energy threshold, 
/ dp s W q u(ps) oc (frw - S2Y 2 - The exponent [i 2 = ^2(9) 
is a function of photon wavenumber and the strength of 
interactions between the bosons. In the relevant limit of 
the weakly interacting gas, the exponent is large [12 S> 1, 
signifying the relative smallness of the photo-solitonic ef- 
fect. As we explain below, such a smallness is associ- 
ated with the quantum orthogonality catastrophe phe- 
nomenon [14| . 

The rest of this paper is organized as follows. In sec- 
tion [IT] we reproduce a derivation of DS solution of the 
Gross-Pitaevskii equation to introduce notations and ter- 
minology. In section IIIII we evaluate the probability to 
excite a specific DS upon absorption of a photon. Section 
HVlis devoted to evaluation and discussion of DSF i.e. the 
total photon absorption rate, resulting in DS formation. 
Finally in section [V] we discuss ways to observe the effect 
experimentally along with the limitations of our theory. 



in Eq. ([7]) and finds two equations for the phase and 
the normalized amplitude x(£)> which are functions of 
£ = x — v s t. The first of these equations acquires the 
form of the continuity relation 



[x 2 



,)]' = o, 



(9) 



where primes denote derivatives with respect to £. Us- 
ing the fact that far from the soliton x(±oo) = 1 and 
i?'(±oo) = 0, one finds i9' = mv s (l — 1/x 2 )- Employ- 
ing this relation, the equation for the amplitude may be 
written in the form 



dU( X ) 



(10) 



where the effective potential U(x), see Fig. [21 is given by 



U(X) 



(i-x 



(11) 



with the Bogoliubov velocity vb = y/cn/r 




(b) 




X 
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FIG. 2: Effective potential U(x) as given by Eq. (TlT) for (a) 
vs > vb; (b) vs < vb- The double-arrowed line designates 
the interval of variation of the normalized amplitude x f° r 
physically allowed solutions of Eq. (|10l) . 



II. DARK SOLITONS 

To establish notations let us briefly discuss the 
localized solutions of the non-linear Gross-Pitaevskii 
equation [ll|. Quasiclassically, this equation is obeyed 
by the condensate wave function: 

id t ^ + ^-d^ + c(n- |*| 2 )* = 0, (7) 
2m 

where n = N/ L is the average concentration and L is the 
length of the system. Hereinafter we switch to the units 
with h = 1. The interaction strength c determines [|[ 
the dimensionless parameter 7 = mc/n whose smallness 
7 <C 1 is the criterion of the weak interaction. 

Looking for a localized solution traveling with a certain 
velocity v s , one substitutes 

*(x,t) = * a (x-v a t) = ^xe i ° (8) 



For v s > vb the potential has the minimum at x — 1 
and the only physically acceptable solutions of Eq. (fTU|) 
are small oscillations around this minimum. In a vicinity 
of x = 1 the potential (|TT|) may be approximated as U ~ 
2m 2 (v 2 — v B )(l — x) 2 an d therefore the small oscillation 
solutions have the form \ — 1 ~ cos{q(x — v s t)) with 
q = 2m\J v"l — v 2 B . Rewriting the last expression as v s — 

^syl + (q/2mvB) 2 , one may recognize it as the phase 
velocity of the Bogoliubov mode. Correspondingly, the 
oscillation frequency qv s = u>b(q) coincides with Eq. ([3]). 
We thus conclude that the only solutions of GP equation 
which travel with a supersonic velocity are Bogoliubov 
quasiparticles. 

The situation is more interesting for v s < vb- In this 
case the potential (jll[) exhibits a maximum at x = 1> a 
minimum at a smaller amplitude and a turning point at 
X = v s /vb < 1. The solution with the proper bound- 
ary conditions, x( = t°°) = 1, is a trajectory which stays 
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at the maximum and then exhibits a bounce down to 
the turning point and back to the maximum. This is 
the DS solution. To find it analytically, one may notice 
that Eq. (fTU)) admits an integral of motion which for DS 
solution reads as 



(12) 



Integrating this equation, one finds for the wave func- 
tion © 



where 



cos i sin — tanh 

2 2 



X — Vet 



cos(6> s /2) = v s /v B 



(13) 



(14) 



with 8 S being the change of phase of the wave function 
across the soliton. The soliton length l s is given by 



i: 



1 = mvB sin(0 s /2) = m^Jv 



(15) 




L X 



FIG. 3: (color online) Density n(x) and phase $(x) profiles 
of a soliton in a system of length L. Note that the density 
perturbation is local, while the perturbation of phase is not. 

The number of particles pushed away from the soliton 
core is 

N s = Jdx(n~\^ s \ 2 ) = ^sin|, (16) 

where the "quantum parameter" K — im/ (juvb) de- 
pends on the inter-particle interaction strength via the 
thermodynamic compressibility which defines the veloc- 
ity vb- Notice that the particle number may be very 
large, N s 3> 1, in the limit of the weakly interacting gas 
(K ^> 1). The energy of the soliton is given by 



dx 



^|^ s | 2 + |(n-|* s |) 2 



4nv B . 3 S 

sin — 

3 2 



(17) 

The calculation of DS momentum requires some care. 
The soliton core momentum, defined by the wave func- 
tion CTJ) is 



p sc — Im dx^ld x ^> = — nsin( 



(18) 



However, one should take into account the periodic 
boundary conditions which ensure that DS phase shift 6 S 



is uniformly spread over the length of the entire system 
L, Fig. [3] Although this does not change the energy of 
the system in the thermodynamic limit (indeed the corre- 
sponding contribution to the energy scales as n9 2 / '(mL)), 
it produces a finite contribution nd s to the momentum. 
As a result the total (core plus the rest of the condensate) 
momentum of the DS state is 



Ps = Psc + 



n9 s = 







(19) 



Equations (flT|) and (fl9|) give an implicit form of DS dis- 
persion relation e s (p s ). The maximum of the soliton en- 
ergy corresponds to 8 S = 7r, where both the soliton ve- 
locity and core momentum vanish v s — p sc = 0. The 
total momentum, however, is finite p s = nn and is uni- 
formly spread across the entire condensate. This is the 
true DS, in a sense that the density vanishes in its cen- 
ter and the particle depletion reaches its maximal value 
N s = 2K/tt. Away from the point 9 S = tt soliton's ve- 
locity is finite v s ^ as well as the density at any point. 
Because of the latter such solitons are sometimes called 
grey. Their velocity approaches sound velocity vb when 
the total momentum approaches zero or 27m, while the 
energy and N s both decrease. Clearly the concept of the 
classical soliton looses sense when the number of particles 
pushed away from the core is comparable to one, N s < 1. 
This takes place when 9 S < 1/K ^ 1, and therefore at 
\Ps\ n/K 3 , and in intervals of the same width around 
the point p = 2im. 

In the limit of the weak interactions 7 <C 1 (i.e. 
K 3> 1) the DS dispersion relation given by Eqs. (fTT)) . 
(fl"9|) . approaches the Lieb II mode, plotted in Fig.[T] The 
convergence is not uniform and the two significantly de- 
viate from each other in the narrow intervals of momenta 
near zero \p s \ < nK~ 3 / 2 and similarly near 27rn [17|. No- 
tice that at the boundaries of this interval the number of 
particles pushed away from the soliton core is still large 
N s tsi \f~K > 1. It is this condition, rather than the 
weaker one N s > 1, which determines the validity of the 
soliton approach. We shall return to this observation in 
section HVl 



III. EXCITATION OF DARK SOLITONS 

Consider a Bose gas subject to a weak space and time 
dependent external potential Vo cos(qx — tot). According 
to the Golden Rule (c.f. Eq. ([2])), the system may absorb 
quanta of this field if its many-body spectrum possesses 
excited states with the momentum q and energy lo. It fol- 
lows from the exactly solvable model Q that such states 
form a continuum whose energy is bound from below by 
the Lieb II mode £2(9)- As argued in the Introduction 
absorption of quanta with the energy in the range given 
by Eq. ([JJ is associated with creation of DS along with 
the phonons or quasiparticles. 

To evaluate the probability of such a process it is conve- 
nient to think of it in terms of the space-time evolution of 
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a state resulting from the photon absorption by the sys- 
tem initially in the ground state. To this end we notice 
that the photon absorption first creates a virtual state of 
the condensate with a local perturbation of the conden- 
sate wavefunction. Since the photon carries momentum 
q and no extra particles, so does the initial local per- 
turbation. Subsequently this perturbation evolves and 
eventually takes a form of a superposition of real excita- 
tions, i.e. conserving overall energy lu in addition to the 
momentum q. We expect that such a final state contains 
a soliton with the momentum p s ~ q and core energy 
s s (p s ) < lu. The small excess energy lu — e s > is carried 
away by phonons, propagating with the sound velocity 
vb- 

The initial separation of the soliton core from a bunch 
of phonons takes a short time, which may be estimated 
as t s = I s /vb- The soliton core is the density depletion, 
which carries momentum — nsmO s (which is very differ- 
ent from p s rs q) and —N s particles. At times t > t s the 
core propagates without dispersion and behaves as a free 
particle with the energy e s . The remaining momentum 
q + ns'm9 s w p s + ns'm6 s = n6 s , cf. Eq. (|19p. and N s 
particles, initially localized on a scale ~ l s , must be car- 
ried away and spread over the entire system at t 3> t s by 
the phonons. As explained above, despite the fact that 
the phonons must carry away large number of particles 
N s and large momentum n6 Sl their final energy is small, 
u> — e s <C VBri6 s . Therefore this is the low-probability 
event, or the "under-barrier" process, which should be 
described as the imaginary time evolution [18|, [l9(] of the 
phonon system (20j . 

To develop such a description we start from the imag- 
inary time r action for the interacting Bose field 



S = I drdx 



2m 



cn\V\ 



-1*1 
2 1 1 



(20) 

It is convenient to parameterize the complex field as 
"J = y/n + (d x ip/n) e lS , where ip(x,r) and $(x,t) are 
two real fields describing density and phase fluctuations 
correspondingly. Assuming small density fluctuations 
d x tp <C Tin and linearizing the resulting action, one finds 



S = Jdr — J dx d x ipd T -d — H s 



(21) 



where the hydrodynamic Hamiltonian of the sound waves 
is given by [2l[ 



VB 

2ir 



dx [K-\d x <pY +K{d x $f 



(22) 



We have omitted terms ~ (d 2 (p) 2 in the Hamiltonian, 
which is equivalent to restricting the spectrum of Bogoli- 
ubov quasiparticles to the phonon branch only. This ap- 
proximation is sufficient for treating photon absorption 
close to the soliton threshold. 

Taking the variations of the imaginary-time action over 



if and $, one finds the semiclassical equations of motion 

d 2 TX ip + iv B Kd 2 J = ttN s [S(x)5{t)-S(x-x)6(t-?)], 
dlJ + i^dly = e s [S(x)5(r)-S(x-x)S(T-f)}. 

The right hand sides of these equations contain sources 
which describe the feedback of DS creation at the point 
x = t = and its subsequent destruction at the point 
x = x, t = f . As discussed above, such creation (destruc- 
tion) of DS is associated with practically instantaneous 
and local injection (removal) of N s particles and momen- 
tum n0 s into (out of) the phonon modes. The sources 
on the r.h.s. of the equations of motion do just that. 
The equations of motions are straightforwardly solved 
by the Fourier transformation. Substituting such a solu- 
tion back into Eq. (|2ip , one finds for the imaginary time 
action 



, , , x — IVbT \ , I „ X + iVBT 

S(x, f) = n + In 1 + tt—— +H- In H 



ih 



where we introduced notations 



iL 



(23) 



(K9 S ± nN s ) 2 K (9 . 9 S \ 2 

^ = a^k = ^ T ± sm T (24) 



and employed Eq. (1 in the last equality in the r.h.s. 
of Eq. (|24|) . Here 9 S is the parameter of a created DS, it 
is related to the DS momentum p s ~q through Eq. ([TO]) . 
The soliton length l s appears in Eq. (|23|) as a short dis- 
tance cutoff. Indeed, one should understand that the ac- 
tual spatial (temporal) extent of the delta-functions on 
the r.h.s. of the equations of motions is the soliton size 
l s (t s =l s /v B )- 

To find a probability of creating DS with the momen- 
tum p s upon absorbing a photon (q,u>), one needs to 
evaluate the Fourier transform of the square of the semi- 
classical matrix element given by e~ s . Specifically, 



W q ,UPs)=Re 



dx dt 

IT 



-S(x,t)—i(q—p s )x-\-i{uj — e s )t 



, (25) 



where we took into account that the momentum p s and 
energy e s (p s ) are carried away by the soliton and there- 
fore should not be absorbed by the phonons. The ana- 
lytical continuation performed in Eq. (|25p to the real fre- 
quencies iu —* lu is accompanied by the time integration 
contour (Wick) rotation r — > it to ensure convergence. 
To evaluate the integral in Eq. (|25"|) , we take into ac- 



count that for small energy excess u> — e s (?) <C mv 2 B the 
range of the allowed soliton momenta p s is rather narrow, 
see Fig. [T] (b), and centered around the photon momen- 
tum q. One may therefore expand the soliton energy as 
£s(Ps) ~ s s (q) + (j) s — q)v s and find for the boundaries of 
the possible soliton momenta, see Eq. (J6|) and Fig. Q] (b) , 



P±(<l,u) = q± 



u) - e s {g) 
vb ± v s (q) 



p+-p-<^ Ps 



(26) 
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Adopting these notations and performing the straightfor- 
ward integrations in Eq. (|25p . one finds 





~P+ - Ps ' 


n+—\ 


Pa - P- ' 




VB 


L Is-' _ 




[ Is' _ 





(27) 



Therefore the soliton creation rate is characterized by the 
power-law dependencies on the deviations of the soliton 
momentum from the upper and lower kinematic bound- 
aries p±(q,io), Fig. [4] (a). The corresponding exponents 
(/i± — 1), see Fig. |4](b), are functions of the soliton param- 
eter 9 S = s (q) and the quantum parameter K as given 
by Eq. (|2"4")) . Since /x_ < fi + , the probability to excite the 
soliton is heavily shifted towards the lower boundary p_ , 
Fig. HJa). he. the soliton is preferentially accompanied 
by the forward moving phonons (in the direction of the 
photon momentum q). 




P- Tiq P+ p s 

T 




FIG. 4: (a) Soliton creation rate as a function of the soliton 
momentum p s . (b) Momentum dependence of the exponents 
H± for K = 10. 



IV. DYNAMIC STRUCTURE FACTOR 

Another quantity of interest is the total absorption rate 
of photons with a given q and lo > £ s {q), which results 
in creation of a soliton with an unspecified momentum. 
This quantity is nothing but DSF S(q, lo) of the ID Bosc 
gas. Integrating W QlU ,(p a ), Eq. (12"?) . over the soliton mo- 
menta p s , one finds for DSF in an immediate vicinity of 
the lower spectral boundary lo > £2(9) 



S(q,io) = I dp s W q ,uip s ) oc — 



P4 



(28) 



where the exponent is given by /12 = fJ-+ + H- — 1. Ac- 
cording to Eq. (|2"4")) . the exponent is expressed through 
the parameters of the soliton 9 S , which in turn is related 
to the soliton momentum through Eq. Q19p , where p s = q. 
As a result 



2K 

M2(g) = — 



sin — 
2 



- 1 



(29) 



For the true DS 8 S = tt and therefore ^{im) sw 0.70JST. 
Notice that at g -> 27m Eq. (J29J) yields /U 2 « 2K - 1, 
different from DSF exponent K — 1 established in the 
framework of the Luttinger liquid theory [2l| . The latter 
is applicable above the dashed line in Fig. Q] 

Employing Eq. (|2"6]) , one may rewrite DSF (|28]) in the 
following form 



S(q,Lo) cc 



1 

't'B 



where the cutoff energy is 



M2(?) 



'(w-e.(?)), (30) 



A(«) = 



—5- sin 3 -i 
2 2 



,(«) 



(31) 



v o 1 



Being multiplied by the intensity of the radiation 
DSF gives a number of solitons excited per unit time 
and per unit length of the irradiated ID gas. 

The power law behavior of DSF near the lower spec- 
tral boundary £2(5) ~ £«(<?) was derived earlier by the 
present authors and M. Pustilnik in Ref. [22J. There a 
mapping between Id Bose and Fermi systems was used 
to prove the presence of the power law non-analyticity 
and evaluate the exponent ^2(9)- However, the method 
adopted there allowed us to deduce the exponent only 
in the limit of strongly interacting bosons 7 » 1 (since 
the latter is mapped onto weakly interacting fermions, 
treated in Ref. [23J ) . Later a method to extract the edge 
exponent \ii for an arbitrary interaction parameter from 
the BA solution was suggested in Refs. [24l. |25|. |26|| . 

Fig. [5] shows comparison between the semiclassical re- 
sult Eq. (f2"9"| and the numerical solution of BA equations 
[23 | for the edge exponent /12 as a function of the mo- 
mentum (in units of n). The agreement between the two 
approaches becomes progressively better for weaker in- 
teractions (the only limit where the soliton picture holds, 
see Introduction). Such an agreement suggests that the 
interpretation of the photon absorption near the lower 
spectral edge as a formation of solitons is indeed consis- 
tent with a fully quantum many-body calculation. The 
latter [2J| does not rely on existence of solitons at all. 
We consider it as a strong confirmation of the thesis that 
absorption of an ac quanta in the frequency window 
results in the formation of DS. 

Notice that even in the weakly interacting limit, the 
semiclassical prediction (f29j) deviates from the exact one 
at very small momenta, see the inset in Fig. [SJa) . This 
is to be expected, since as was discussed in section [TTl 
the soliton picture looses its validity at sufficiently small 
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FIG. 5: (color online) Edge exponent ^(q), of. Eq. (|30p . as a 
function of momentum in units of n for (a) 7 = 0.05, K = 33; 
(b) 7 = 0.4, K = 5.2; (c) 7 = 1, K = 2.7. The full (blue) line 
is the semiclassical result (|29[) ; dashed (red) line is the Bethe 
ansatz solution of Ref. 12411. 



momenta. Inspecting Eq. (|29[) , one notices that the semi- 
classical exponent becomes negative at 9 S < n/y/K, con- 
trary to the exact results, Fig. [5l Using Eq. (fl6|) . one 
finds that this corresponds to the number of missing par- 
ticles in the soliton core N s w y/K. This observation col- 
laborates with the discussion presented in the end of sec- 
tion mi which suggests that the semiclassical treatment 
looses validity for the very grey solitons with N s < \f~K, 
i.e. \p a \ < nK~ 3 / 2 , cf. Eq. (HI]). We stress that the 



power-law behavior of DSF at the exact lower spectral 
boundary £2(9) is valid for any momentum. However for 
q < nK~ z / 2 the semiclassical approximation for the ex- 
ponent [ii (q) fails. Instead, the exact exponent [13, l2~j ] 
scales linearly with momentum 112(1) ~ K 3 / 2 q/n. 



V. DISCUSSION 

We have shown that the absorption of a photon with 
the energy above a certain threshold leads to formation 
of the DS. In the narrow energy window e s < uj < e s + 
A the total soliton formation rate per unit length of ID 
Bose cloud is given by S(q,ui)VQ, where DSF is given by 
Eq. ((30]) and Vq cos(ujt — qt) is the external ac potential 
applied to ID Bose gas. The momentum-resolved rate is 
given by Eq. (|27p . i.e., solitons with the momentum in 
the interval p s ± dp s /2, where p s belongs to the window 
©, are created with the rate W q , u (p a )Vodp a . 

An important question is what are the corresponding 
rates for a larger energy of the photon: e s + A < lu < vsq. 
According to the arguments given in the Introduction, 
absorption of such a photon should necessarily lead to 
DS formation. Yet our calculations are not directly ap- 
plicable in this case. Indeed, we have used linearized 
dispersion relation for the quasiparticles (phonons) ex- 
cited along with DS. For energies above e s + A such an 
approximation is not valid. This is because Bogoliubov 
quasiparticles with momenta above mvs that take the 
excess energy can not be well approximated by phonons. 
The photon absorption is dominated by creation of a DS 
and single quasiparticle moving in the direction of the 
wavevector q (i.e. moving forward). Thus the parameters 
of the typical DS may be found by plotting a replica of 
the Bogoliubov spectra which starts at some point along 
the absorption edge e — e s (q) and passes through (q,w). 
The starting point prescribes DS momentum and energy. 
We expect that the power-law Eq. (|30)) saturates at an 
excess energy of order A(q), i.e. at uj — e 2 ~ A. As a re- 
sult, DS formation rate per unit length may be estimated 
as (Vq /h 2 VB)e~ atJ ' 2 ^ q \ with a numerical factor a ~ 1. 

Recalling that for a true DS ^ 2 = 0.7AT = l.liV a , 
one realizes that the photo-solitonic rate is exponen- 
tially suppressed with the increase of DS depleted parti- 
cle number N s . Physically the origin of this smallness is 
in the orthogonality phenomenon: the state of the sys- 
tem immediately after absorption of the photon is al- 
most orthogonal to the state with the soliton causing a 
re-distribution of density and phase of the condensate 
in the one-dimensional system. The corresponding ma- 
trix element is exponentially small in the parameter N s . 
This fact dictates a rather stringent limitations on the 
experimental observability of the photo-solitonic effect. 
Increasing interactions (i.e. decreasing K) makes the ex- 
ponential factor in the photo-solitonic rate less severe, on 
the other hand it simultaneously decreases N s , making 
it more difficult to observe the excited solitons. Assum- 
ing Vo/[2nK) = 100 Hz, Ref. Q we estimate the soliton 
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production rate for the relatively "light" solitons with discussions. This research is supported by DOE Grant 
N s = 10 as 1 event per Bragg pulse of duration of 2 No. DE-FG02-08ER46482 and A.P. Sloan foundation, 
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Acknowledgments 

We thank A. Abanov, J.-S. Caux, D. Gangardt, D. 
Gutman, V. Gurarie, and A. Imambekov for numerous 



[1] J. Denschlag, et al, Science 287, 97 (2000); S. Burger et 

al., Phys. Rev. Lett. 83, 5198 (1999); B.P. Anderson et 

al, Phys. Rev. Lett. 86, 2926 (2001). 
[2] R. Dum, et al, Phys. Rev. Lett. 80, 2972 (1998). Y. S. 

Kivshar, B. Luther-Davies, Phys. Rep. 298, 81 (1998). 

A. Muryshev et al, Phys. Rev. Lett. 89, 110401 (2002). 
[3] S. Burger et. al, Phys. Rev. Lett. 83, 5198 (1999). B. 

Wu, J. Liu, and Q. Niu, Phys. Rev. Lett. 88, 034101 

(2002). 

[4] P. Lenard, Uber die hchtelektrische Wirkung, Annalen 

der Physik, 8 (1902). 
[5] A. Einstein, On a Heuristic Viewpoint Concerning the 

Production and Transformation of Light, Annalen der 

Physik, 17 (1905). 
[6] J. Stenger et al, Phys. Rev. Lett. 82, 4569 (1999). 
[7] J. Steinhauer, R. Ozeri, N. Katz, and N. Davidson, Phys. 

Rev. Lett. 88, 120407 (2002); J. Steinhauer et. al., Phys. 

Rev. Lett. 90, 060404 (2003) and references therein. 
[8] In Bragg scattering tecnique strictly speaking, we talk 

about two-photon absorption. 
[9] E.H. Lieb and W. Liniger, Phys. Rev. 130, 1605 (1963); 

E.H. Lieb, Phys. Rev. 130, 1616 (1963). 
[10] N. A. Slavnov, Teor. Mat. Fiz. 79, 232 (1989); Teor. Mat. 

Fiz. 82, 389 (1990). 
[11] J.-S. Caux and P. Calabrese Phys. Rev. A 74, 031605 

(2006). 

[12] E.M. Lifshitz and L.P. Pitaevskii, Statistical Physics, 

Part 2 (Pergamon Press, 1980). 
[13] P. P. Kulish, S. V. Manakov, L. D. Faddeev Theor. Mat. 

Fiz. 28, 38 1976. 



[14] P. W. Anderson, Phys. Rev. Lett. 18, 1049 (1967). 

[15] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, 

Rev. Mod. Phys. 71, 463 (1999). 
[16] T. Tsuzuki, J. Low Temp. Phys. 4, 441 1970. 
[17] V. Gurarie, D. Gangardt, M. Khodas, A. Kamenev, and 

L.I. Glazman, unpublished. 
[18] L.D. Landau and E.M. Lifshitz, Quantum Mechanics, 

Non-Relativistic Theory (Pergamon Press, 1977). 
[19] L.S. Levitov and A.V. Shytov, Pis'ma Zh. Eksp. Teor. 

Fiz. 66, 200 (1997) [JETP Lett. 66, 214 (1997)] . 
[20] S. Brazovskii and S. I. Matveenko, Sov. Phys. JETP 96, 

555 (2003) [Zh. Eksp. Teor. Fiz. 123, 625 (2003)]; S. 

Brazovskii and S. I. Matveenko, Phys. Rev. B 77, 155432 

(2008). 

[21] V.N. Popov, Theor. Math. Phys. 11, 565 (1972); K.B. 
Efetov and LA. Larkin, Sov. Phys. JETP 42, 390 (1975) 
[Zh. Eksp. Teor. Fiz. 69, 764 (1975)]; F.D.M. Haldane, 
Phys. Rev. Lett. 47, 1840 (1981). 

[22] M. Khodas, M. Pustilnik, A. Kamenev, and L.I. Glaz- 
man, Phys. Rev. Lett. 99, 110405 (2007). 

[23] M. Pustilnik, M. Khodas, A. Kamenev, and L. I. Glaz- 
man, Phys. Rev. Lett. 96, 196405 (2006). 

[24] A. Imambekov and L.I. Glazman, Phys. Rev. Lett. 100, 
206805 (2008). 

[25] R. G. Pereira, S. R. White, and I. Affleck, Phys. Rev. 
Lett. 100, 027206 ( 2008); V. V. Cheianov and M. Pustil- 
nik Phys. Rev. Lett. 100, 126403 (2008). 

[26] M. Khodas and A. Kamenev, unpublished. 



